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Accounting  for  the  Penetration  Entrance  Effect 

Yehuda  Partom 


Abstract 

By  running  simulations  in  which  we  constrained  the  backward  motion  of  the 
entrance  boundary,  we  show  that  the  entrance  effect  is  not  caused  by  the  boundary 
motion.  We  use  a  simple  time  dependent  Tate  model  to  show  that  the  entrance  effect 
is  caused  by  the  initial  (ID  strain)  high  penetration  velocity,  and  by  the  time  depen¬ 
dency  of  the  penetration  process. 

1.  Introduction 

Running  a  computer  simulation  of  a  long  rod  penetration  event,  one  can  see 
clearly  an  entrance  boundary  effect.  Figures  1  and  2  are  output  plots  from  an  AUTO- 
DYN/EULER  run  of  a  L/D  =  1 0  tungsten  alloy  rod  at  V  =  1 .5  km/s  into  a  RHA  steel  tar¬ 
get.  The  figures  are  taken  from  our  report  about  the  projectile-flow  effect  [1].  The 
plots  in  Figures  1  and  2  are  from  a  run  with  a  strain  softening  projectile  material 
according  to: 

Y=Y0(i-pey  ,  (1) 

where:  Y  =  flow  stress,  Y0  =  yield  stress  =  2  GPa,  e§ff  =  effective  plastic  strain,  and 
(3  =  0.25.  We  chose  to  show  the  results  of  this  run  because  strain  softening  enhances 
the  contrast  between  the  entrance  and  the  quasi  steady  state  regimes. 

We  see  that  initially,  penetration  velocity  is  significantly  higher  than  during  the 
quasi  steady  state  phase.  Also,  deceleration  of  the  penetration  velocity  is  higher  up  to 
t  =  25  (is  (penetration  of  about  three  diameters).  This  behavior  is  typical  of  long  rod 
penetration  events,  and  we  refer  to  it  as  the  "entrance  effect."  In  contrast  to  what  we 
get  from  a  simulation,  there  is  no  entrance  effect  in  a  Tate  model  run. 

In  Figure  3,  we  show  a  penetration-erosion  curve  from  a  Tate  run  of  the  same 
problem,  with  Yp  =  2  GPa  and  R(  =  5  GPa.  We  see  no  entrance  effect,  and  the  end 
(secondary)  phase  also  looks  different. 
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TUNGSTEN  INTO  STEEL 


Figure  1.  Penetration-erosion  curve  from  an  AUTODYN  run  of  a  L/D  =  10  tungsten  alloy 
rod  penetrating  RHA  steel  at  1.5  km/s. 


TUNGSTEN  INTO  STEEL 

Figure  2.  Penetration  velocity  and  projectile  tail  velocity  histories  for  the  same  run  as  in 
Figure  1 . 
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Tate  run 


Figure  3.  Penetration-erosion  curve  from  a  Tate  run  of  a  L  =  100  mm  tungsten  alloy  rod 
into  RHA  steel  at  1 .5  km/s.  Yp  =  2  GPa,  Rt  =  5  GPa 


In  Figure  4,  we  show  a  typical  plot  from  an  AUTODYN  simulation  of  material 
boundaries  and  velocity  arrows  at  about  9  jlls  (~  YiD)  into  the  penetration  process.  We 
see  that  a  typical  lip  is  formed  on  the  target  front  surface  around  the  crater.  Because  of 
this  lip,  and  the  backward  target  motion  that  creates  it,  we  used  to  think  that  target 
resistance  is  lower  near  the  entrance  boundary,  leading  to  the  entrance  effect.  We 
used  to  advocate  this  "understanding"  quite  vigorously,  and  even  devised  a  way  to 
estimate  the  reduced  resistance  using  a  cavity  expansion  approach  [2],  This  view  is 
not  shared  by  other  investigators  [3]  who  claim  that  the  entrance  effect  results  from  the 
initial  shock  phase  and  the  high  (1 D  strain)  particle  velocity  associated  with  it. 
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Figure  4.  Material  boundaries  and  velocity  arrows  plot  from  an  AUTODYN  simulation  at 
about  Vi D  penetration. 


In  what  follows,  we  investigate  the  origin  of  the  entrance  effect  by  means  of 
AUTODYN  simulations.  We  do  this  by  inhibiting  the  backward  motion  of  a  part  of  the 
entrance  boundary.  We  find  that  the  entrance  effect  is  not  affected  by  the  inhibition  of 
the  backward  motion  of  the  entrance  boundary. 

Recently  Walker  and  Anderson  [4]  proposed  a  time-dependent  model  for  long 
rod  penetration,  so-called  WA  model.  It  is  based  on  integration  of  the  momentum 
equation  along  the  projectile-target  center  line.  Among  other  things,  Walker  and 
Anderson  show  that  their  model  captures  quite  well  the  entrance  effect. 

In  what  follows,  we  use  Walker  and  Anderson's  approach  to  formulate  what 
can  be  called  a  time  dependent  Tate  model.  It  is  a  much  simplified  version  of  the  WA 
model.  Running  this  model,  we  show  that  it  too  is  able  to  capture  the  entrance  effect. 
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2. 


Simulations 


We  are  using  AUTODYN/EULER  version  2.65.  The  tungsten  alloy  (projectile) 
and  RHA  steel  (target)  material  models  and  parameters  are  the  same  as  in  [1]  and  [5]. 
The  projectile  radius  is  D/2  =  5  mm,  and  the  Euler  cell  size  is  1  x  1  mm  so  that  there 
are  five  cells  across  the  radius.  This  may  not  be  enough  to  obtain  convergent  results, 
as  shown  in  [6].  But  as  our  purpose  is  to  investigate  the  origin  of  the  entrance  effect 
(and  not  to  compare  simulations  to  experiments),  five  cells  across  the  radius  seem  sat¬ 
isfactory.  The  projectile  aspect  ratio  is  L/D  =  1 0,  and  the  impact  velocity  V  =  1 .5  km/s. 

We  terminated  our  runs  after  400  cycles  (approximately  40  ps)  and  plotted 
material  status  plots  every  100  cycles. 

In  all  runs  the  entrance  boundary  is  stress  free  on  a  ring  D/2  <  r  <  D/2  +  Ar,  and 
constrained  in  the  axial  direction  elsewhere  (r  >  D/2  +  Ar).  The  different  runs  are  list¬ 
ed  in  Table  1 . 


Table  1 


List  of  Simulation  Runs 


No. 

Ar  mm 

1 

i 

2 

2 

3 

5 

4 

10 

5 

oo 

Results  for  a  regular  run  (unconstrained  boundary,  run  no.  5)  are  shown  in 
Figures  5  and  6).  In  Figure  5,  we  show  material  status  plots  at  100,  200,  300,  and  400 
cycles.  In  Figure  6,  we  show  velocity  vectors  near  the  interface  at  the  same  times. 

Figure  7  shows  material  status  plots  at  cycle  300  for  runs  1  to  4,  and  in 
Figure  8,  the  corresponding  velocity  vectors  plots.  We  see  that  when  Ar  =  1  and  2 
mm,  the  target  material  is  not  able  to  flow  through  the  gap.  For  Ar  =  5  and  10  mm, 
target  material  flows  through  the  gap  and  a  lip  is  formed.  But  in  all  cases,  the  penetra¬ 
tion  is  the  same,  and  the  mushroom  configuration  is  the  same.  The  amount  of  back¬ 
wards  flow  of  the  target  boundary  does  not  seem  to  influence  the  penetration  process. 
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Figure  5.  Material  status  plots  for  a  regular  run  (no.  5). 
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Figure  6.  Velocity  vectors  plots  for  a 
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Figure  7.  Material  status  plots  at  300  cycles  for  runs  1  to  4. 
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Figure  8.  Velocity  vectors  plots  at  300  cycles  for  runs  t  to  4. 
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To  see  this  even  more  clearly,  we  show  penetration  velocity  histories  and  pene¬ 
tration-erosion  curves  for  runs  1  and  5  in  Figures  9  and  10  .  From  Figure  9,  we  can 
see  that  the  penetration  velocity  histories  are  practically  the  same.  From  Figure  1 0,  we 
can  see  that  the  penetration-erosion  curves  are  also  practically  the  same,  and  that  they 
both  show  the  entrance  effect. 

We  conclude  from  the  simulations  that  the  entrance  effect  is  not  caused  by  low 
target  resistance  near  the  entrance  boundary. 

In  the  next  section  we  propose,  as  did  Walker  and  Anderson  [4],  that  the 
entrance  effect  can  be  accounted  for  by  time  dependency  of  the  penetration  process. 

3.  Time  Dependent  Tate  Model 

Walker  and  Anderson  [4]  have  shown  that  by  integrating  the  axial  momentum 
equation  along  the  projectile-target  center  line,  one  gets  the  following  relation: 

Ppf  ft^X  +  Ptf  ^7dx  +  2  pp(Cu(xi)]2-[u(xp)]2}  +  ypt  ([u(xt)]2-  [lj(x  j)]2  ) 

/xn  Jx-,  Z 


~  [°XX  "  °xx  (*p)  ]  "  2 

where: 

x  =  axial  direction,  y  =  transverse  direction,  u  =  axial  particle  velocity, 
t  =  time,  Ojj  =  cartezian  stress  components,  pp  =  projectile  density  (assumed 
constant),  pt  =  target  density  (assumed  constant),  xp  <  x  <  xt  is  the  integration 
range  along  the  center  line  (y  =  0). 

Unlike  Walker  and  Anderson,  we  choose  (in  what  follows)  xp  to  be  the  elastic- 
plastic  boundary  in  the  projectile,  and  xtthe  elastic-plastic  boundary  in  the  target.  Xj 
is  the  projectile-target  interface. 


9a 


xy 


9y 


dx 


dx  =  0 


(2) 
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In  what  follows,  we  use: 


u(Xp)  =V  (projectile  tail  velocity) , 

(3) 

u(Xj)  =  U  (penetration  velocity) , 

(4) 

and  assume: 

C 

III 

o 

(5) 

Walker  and  Anderson  went  on  to  evaluate  the  last  two  integrals  in  Eq.  (2)  by 
assuming  a  specific  form  for  the  flow  field  u(x,  y).  We  are  using  a  simpler  approach. 
Recognizing,  as  did  Walker  and  Anderson,  that: 

Pdaxy 

'°“(X,)'2  ~3fdx  =  R' 

Jx  i 

(6) 

P  aaxy 

axx^V  2  1  3y  C*X-YP' 

Jxp 

(7) 

where  Rt  and  Yp  are  the  Tate  model  target  resistance  and  projectih 
tively. 

s  strength,  respec- 

We  get  from  Eq.  (2): 

Pt(  ffd*'2  P,u2  +  Rt  =  -Pp(  '^dx-lppl^-vVp 

Jx\  Jx  p 

(8) 

To  evaluate  the  integrals  in  Eq.  (8),  we  make  the  following  simplifying  assump¬ 
tions: 

Xj  -  xp  =  const.  =  Sp  , 

(9) 
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xt  -  Xj  =  const.  =  St , 


(10) 


u  =  U  +  (V-U)  fp  (x)  in  the  projectile,  (1 1 ) 

and 

u  =  Uft  (x)  in  the  target,  (1 2) 

where  the  functions  f  and  ft  are  advected  without  change  at  a  velocity  U,  as  shown  in 
Figure  11 . 


u 


Figure  1 1 .  Simplified  flow  field  on  center  line. 

The  advection  assumption  leads  to: 

dfp  _  ydfp 

8t  dx 


(13) 
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aft=.u3fj 

3t  3x 


so  that  in  the  projectile: 


fudu+(dy.du)f U(V.U) 

3t  dt  v  dt  dt  /  p 


afp 

3x 


and  in  the  target: 


5u  _  dU  f  .  .2  3ft 
¥'  dt  ,_u  ¥ 


Evaluating  the  integrals  in  Eq.  (8)  we  get: 


^dx  =  K,s,dU+u2 

at  dt 


u(v-u) 


where  Kp/  iq  are  shape  factors.  For  fp/  ft  linear,  Kp  =  Kt  =  1/2. 
degree  parabolas,  Kp  =  iq  =  1/3. 

Substituting  for  the  integrals  from  Eqs.  (17)  and  (18)  into 
ing,  we  finally  get: 


Rt  +  “Pt  U2  +  [Pp  Sp  (1  -  Kp  )  +  pt  St  Kt]  - 
=  Yp  + 1-  pp  (V-U)  2  -  pp  Sp  Kp  ^ 


(14) 


(15) 


(16) 


(17) 


(18) 

For  fp,  ft  concave  2nd 
Eq.  (8)  and  rearrang- 


(19) 
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We  see  that  Eq.  (19)  collapses  to  Tate's  equation  (modified  Bernouli  equation)  when 


dt^  _<dt  (steady-state  assumption).  We  therefore  refer  to  Eq.  (19) 

as  the  "time 

dependent  Tate  equation."  Using  Eq.  (19)  together  with  the  kinematic  relations  of 

Tate's  model  we  thus  define  a  time  dependent  Tate  model,  as  follows: 

dU=F(u,V,^) 

(20) 

as  defined  by  Eq.  (19). 

dV  _  Yp 

dt  Pp  (L  -Sp) 

(21) 

4^  =  -(V-U) 
dt 

(22) 

dp  =  u 

dt 

(23) 

with  the  initial  conditions: 

V  (0)  =  VQ  (impact  velocity) 

(24) 

U  (0)  =  U0  (ID  strain  interface  velocity) 

(25) 

L  (0)  =  Lc 

(26) 

p  (0)  =  0 

(27) 

We  integrate  the  four  ODEs  Eqs.  (20)  to  (23)  using  a  standard  ODE  system  solver. 

At  the  end  phase  of  the  integration,  |  dV/dt  |  becomes  large  as  L  -  Sp  becomes 
small.  Because  of  that,  using  Eq.  (19)  to  evaluate  dU/dt  becomes  unreliable.  We 
therefore  switch  to  a  more  robust  way  to  compute  dU/dt  as  follows: 

/dU) 

We  compute  \  ^  f_^  by: 
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(28) 


fdU)  =  pP  (V'U)  dV 
l  dt  lJate  ptU  +  pp  (V-U)  dt 

as  obtained  from  differentiating  Tate's  equation  and  we  get: 


(-dU)  =  max 

(-dU) 

|.dU) 

dt  ' 

dt  Eq. 1 9 

dt  Tate] 

To  see  how  the  time  dependent  Tate  model  performs,  we  compare  it  to  the  reg¬ 
ular  Tate  model  and  to  simulation  results.  The  example  we  show  is  for  an  L/D  =  10 
tungsten  alloy  rod  penetrating  an  RHA  steel  target.  The  model  parameters  are: 

Pp  =  1 7-3  g/cc,  Yp  =  2  GPa,  Sp  =  1 4mm,  Kp  =  0.5, 

pt  =  7.85  g/cc,  Rt  =  6  GPa,  S{  =  65  mm,  Kt  =  0.3, 

VQ  =  1 .5  km/s,  U0  =  1  km/s,  L0  =  1 00  mm. 

In  Figures  12.1  to  12.3,  we  show  history  plots  from  the  simulation  (U(t),  p(t), 
and  p(L)). 

In  Figures  13.1  to  13.3,  we  show  the  same  history  plots  obtained  from  the  time 
dependent  Tate  and  from  the  regular  Tate  models.  We  copied  on  these  plots  the 
curves  from  Figures  12.1  to  12.3,  respectively. 

We  see  the  following: 

-  The  entrance  effect  is  captured  quite  well  by  the  time  dependent 
Tate  model. 

-  As  expected,  the  time  dependent  Tate  does  not  do  well  at  the  end  phase. 
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Figure  1 2.1 .  Penetration  velocity  history  from  the  simulation. 
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Figure  12.2.  Penetration  history  from  the  simulation. 
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Figure  12.3.  Penetration-erosion  history  from  the  simulation. 
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Figure  1 3.1 .  Penetration  velocity  histories. 
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Figure  13.2.  Penetration  histories. 


penetration  erosion 


Figure  13.3.  Penetration-erosion  histories. 
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4. 


Conclusions 


We  ran  computer  simulations  with  AUTODYN  to  investigate  the  origin  of  the 
entrance  effect  in  long  rod  penetration.  By  constraining  the  backward  motion  of  the 
entrance  boundary,  we  show  that  the  entrance  effect  does  not  depend  on  the  motion 
of  the  entrance  boundary. 

Based  on  Walker  and  Anderson's  approach  [4]  of  integrating  the  momentum 
equation  along  the  axis,  we  propose  a  time  dependent  Tate  model.  Running  the 
model  and  comparing  to  simulation  results  we  show  that  the  model  captures  quite 
well  the  entrance  effect. 

We  conclude  that  the  entrance  effect  results  from  the  initial  high  (1  D  strain) 
penetration  velocity,  and  from  the  time  dependency  of  the  penetration  process. 
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